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. We investigate invariants for random elements of different hyperbolic 
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groups. We provide a method, using Cayley graphs of groups, to compute 



Q ' the probability distribution of the minimal length of a random word, and ex- 

plicitly compute the drift in different cases, including the braid group -B3. We 
also compute in this case the return probability. The action of these groups 
on the hyperbolic plane is investigated, and the distribution of a geometric 
invariant, the hyperbolic distance, is given. These two invariants are shown 

> ■ 

■ to be related by a closed formula. 
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O I. INTRODUCTION 

o 

Q-r This paper is devoted to a systematic study of random walks on the modular group 

PSL{2, Z) and on some closely related groups: the braid group B 3 , the Hecke groups H q 
and the free groups F n (all definitions are given below). We study simultaneously the limiting 
distribution of random walks on the Cayley graphs of these groups and on their Riemann 
surfaces. We analyze the statistical properties of random walks on the Cayley graphs of 
the above mentioned groups both in a metric of words and in the natural metric of the 
hyperbolic plane. 

The very subject of our investigation is not new — the statistics of Markov chains on the 
subgroups of the group PSL(2, IR) has been extensively studied in the mathematical liter- 
ature. Among the known results connected to the theme of our work we can mention: (a) 
the central limit theorem for Markov multiplicative processes on discrete subgroups of the 
group PSL(2, IR) jn],|I7]], (b) some particular examples of exact results for limiting distribu- 



tion functions of random walks on Cayley graphs of free and modular groups |11|J2^J25|J22 
and (c) conjectures concerning the return probability and drift on the braid group B 3 

In the present work we rigorously compute the drift and the return probability for sym- 
metric random walks (in metric of words) on the groups H q and B 3 . Moreover, as it has 
been said, we pay a special attention to the statistics of random walks on the Riemann sur- 
faces of the groups PSL(2, Z), B 3 , H q , F n . Namely, we study a 2 x 2 matrix representation 
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of these groups and consider their homographic actionP] on the hyperbolic plane 7i. This 
allows to embed the Cayley graphs in Ti, and to define isometric hyperbolic lattices. Taking 
advantage of the hyperbolic metric on Ti, we investigate the probability distribution of the 
geodesic distance between ends of random processes with symmetric transition probabilities 
on these lattices of 7i. We show that this problem reduces to the study of the modulus of a 
random product of matrices. The part of our investigation is semi-analytic and is based on 
numerical results on the structure of the invariant distribution of geodesies at the boundary 
of Ti. We found very interesting the fact that the drift on a Cayley graph in a metric of 
words coincides after proper normalization with the drift on the corresponding isometric 
lattice of Ti in the natural hyperbolic metric. This result establishes a nontrivial relation 
between two group invariants: in one hand the irreducible length of an element, which does 
not depend on the representation, and on the other hand, the hyperbolic distance associated 
to an element (directly linked to its modulus), defined only for this matrix representation. 

As an application of our results, we consider the relation between the distribution of 
Alexander knot invariants and the asymptotic behavior of random walks over the elements 
of the simplest nontrivial braid group B 3 . This class of problems arises naturally even 
beyond the aims of our particular investigation: the limiting behavior of Markov chains on 
braid and so-called "local" groups can be regarded as a first step in a consistent development 
of harmonic analysis on branched manifolds (Teichmiiller spaces are an example). 

The paper is structured as follows. In section |I| we give the basic definitions and in- 
troduce the different groups and their Cayley graphs. A general solution of the diffusion 
problem on these graphs, as well as exact computations of the drift and the return proba- 



bility for B 3 are developed in section |TT[ Section [IV| is devoted to the study of the action 
of these groups in the hyperbolic plane; a discussion of our results and the relation between 
the different approaches are presented in section [V|. 

II. HYPERBOLIC GROUPS PSL{2,Z), H Q , B 3 , F N AND THEIR CAYLEY 

GRAPHS 

A. Basic definitions 

We consider a special class of so-called hyperbolic groups - the modular group PSL(2, Z), 
and some of its generalizations - the Hecke groups H q and the braid group B 3 . We also 
recall already known properties of the free groups F n , usefull in the context of our work. 



1 since the 2-representation of B 3 is not unimodular, this action is not faithfull 
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1. The modular group PSL(2, Z) is a free product Z 2 * Z 3 of two cyclic groups of 2nd 
(generated by a 2 ) and 3rd (generated by 63 ) orders. In a standard framing using generators 
S (inversion) and T (translation), the group PSL(2, Z) is defined by the following relations 

(srr = bi = i m 

S 2 = a\ = 1 1 ] 

Being a discrete subgroup of the group PSL(2, IR), the generators T and S of the modular 
group PSL(2, Z) have a natural representation by unimodular matrices T and S: 



r=LJi , J ( 2 ) 





2. In addition to the modular group PSL(2, Z) we shall consider the so-called Hecke 
group H q which "interpolates" between the modular group (for q = 3) and the free group 
F3 with 3 generators, the so-called A group (for q = 00). The Hecke group H q is isomorphic 
to Z 2 * Z 9 (we denote by a 2 and 6 9 the generators of orders 2 and g). It is defined by 
straightforward generalization of the relations (|T]) 

(ST q )« = 6g = l , . 

5 2 = al = 1 1 j 

and the generators T q and 5 have the following matrix representation (compare to (0)): 

The parameter q takes the discrete values q = 3, 4, 5, 6, . . .. 

3. The braid group B 3 is defined by the following commutation relations among genera- 
tors {a h a 2 }: 

a 1 a 2 cr 1 = cr 2 criO- 2 /c x 

-1 -1 ( 5 ) 

(Ji(Ji = cr 2 cr 2 = e 

In our further construction we shall repeatedly use the following framing: 

r -1 -1 (6) 
b = a 1 a 2 

The generators of the group B3 can be represented by PGL(2, IR)-matrices. To be more 
specific, the generators o\ and <r 2 in the Magnus representation || read 
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where t is the free parameter. We conveniently introduce the parameter u = y— t, and 
consider normalized generators of determinant 1: 



»^ ; «* 3 = A/ttU (8) 





The group generated by u <j\ and M a 2 will be denoted later on as PSL(2, Z) u . Indeed it is 
just a deformation of PSL(2, Z), which preserves all its commutation relations. For u — 1, 
one has PSL(2, Z) u = PSL(2, Z), and the group _B 3 is a central extension of PSL(2, Z) u of 
center 

(axa 2 ) 3A = (^) 3A = (a^ 2 ^) 2A = (a 2 a 1( r 2 ) 2A = f ^ J , VA E Z (9) 

(let us note that the center is isomorphic to Z) . Recall that graphically, to each word of B 3 
correspond a particular three-strand braid, going from above downwards, (see Fig.|l]). A 
closed braid is obtained by gluing the "top" and "bottom" free ends on a cylinder. Any closed 
braid defines a link (in particular, a knot). However the correspondence between braids and 
knots (links) is not one-to-one and each link (knot) can be represented by infinite number of 
different braids (see The irreducible length of a braid gives nevertheless an interesting 

characteristic of the link complexity. 
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FIG. 1. Some examples of closed braids and corresponding links. 
There exists an extensive literature on general properties of braid groups — see [H; for 



the last works on the normal forms of words, we shall quote [10]. 

Any element of the group G = {PSL(2, Z), H q , B 3 } is defined by a word in alphabets of 
corresponding letters (generators): 

• {S, T, T -1 } or {a 2 ,b 3 ,b^}— for PSL(2,Z) 

• {S^q^T- 1 } or {a 2 ,b q ,b- 1 }— for H q 

• {cr 1 ,a 2 ,(Ti 1 ,(T2 1 } or {a, a" 1 , 6, for B 3 . 

We denote by w n a word corresponding to a given record of length n, and by L (w n ) the 
irreducible length in the metric of words (the superscript G is precised only when it is 
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necessary), or in other terms the minimal number of generators necessary to build w n . The 
irreducible length can be also viewed as a distance from the unity on the Cayley graph of 
the group G. Note that L G {w) depends on the set of generators we consider. 

B. Cayley graphs 

The modular group PSL(2, Z) is a particular case of the Hecke group H q at q = 3. We 
therefore consider without any loss of generality the Cayley graphs of the groups H q for 
q = 3,4, .... The Cayley graph of B3 will be constructed afterwards. We investigate in this 
part only the abstract presentation of the groups in terms of commutation relations and do 
not pay attention to any representation. We recall that the Cayley graph of a group G is the 
graph whose vertices are labeled by group elements, and whose links are as follows: w and 
w' are linked if and only if there exists a generator g such that w' = wg. Following this rule, 
we can easily construct the Cayley graph Q q of the group H q represented by {a,2,b q ,b q 1 }. 
For any finite values of q the graph Q q has local g-cycles (because b q is of order q), while 
the corresponding dual (or "backbone") graph is the tree graph ~W q , which is precisely the 
graph of F q . This is due to the free product structure of H q ~ Z 2 * Z g (see explanations 
below). The graph Q q is shown in Fig.|2], where the backbone graph is marked by a dotted 
line. 




FIG. 2. Cayley graph of H q , with q = 4. Arrows correspond to generator b q = ST q , thin lines 
to generator 02 = S = S~ x . The dashed graph is the corresponding backbone graph, a tree T g , 
graph of F q . The distance d between A and B is d = 5, while the distance k (corresponding to the 
distance between A' and B') is k = 2. 

III. DIFFUSION ON GRAPHS 

In this section we investigate some statistical properties of random walks on the groups 
introduced above, using their Cayley graphs. In particular we consider simple random walks, 
that are walks of nearest neighbour type with symmetric transition probabilities. 
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A. Random walk on PSL(2, Z) and H q 



We consider free product groups of the form Z2 * Z g (isomorphic to H q ), in the framing 
which uses the generators 02 and b q , of order 2 and q respectively. Graphs of such groups 
are shown in Fig.|2]. We define two different "metrics" of words on those graphs. The first 
metric is associated with the geodesic distance d on the graph — the minimal number of steps 
between two points, and the second metric is associated with the geodesic distance k (called 
later the "generation") on the backbone graph ~W q . Our goal is to compute the probability 
P q (d, n) of being at a distance d from the initial point (the root of the graph) after n random 
steps. The probability P q (k,n) of being on the backbone graph at a generation k from a 
root point after n random steps will also be of use. 

First of all we compute P 3 (d,n) and P 3 (k,n) for the case of PSL(2, Z). In this case the 
graph structure ensures the relation 

P 3 (2x,n) = P 3 (x,n) 

Therefore we can consider only P 3 (k,n). Write 

P 3 (k,n) = P 3 i (k,n) + Pg{k,n) 

distinguishing for an elementary triangular cell located at generation k the vertex closest 
to the root (corresponding to P 3 (k,n)) and the two others (corresponding to P 3 (k,n)) (see 
Fig.|7|). A direct enumeration gives the following master equation for k > 2: 



n + l) = \ (P°(k, n) + P°(k - 1, n 
P 3 °(k, n + l) = lP*(k, n) + § (P£(k, n) + P&k + 1, 



(10) 



with initial conditions of the form 

P 3 l (k,0)=a5 kfi (n) 
P 3 °(M) = (l-a)5 M 

where a is an arbitrary parameter fixing the initial condition and varying in the interval 
< a < 1. We are seeking for the asymptotic (1 <C fc < n) solution to fliTf ) near the 
maximum of the probability distribution, and therefore will not take into account the specific 
form of the boundary condition. 

Define the Laplace-Fourier transform: 

00 +00 

Q^ix, s)=T [Pr] = s n £ e ikx P 3 '°(k, n) (12) 

n=0 k=— 00 
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whose inverse can be written in the form 

1 r ds 

F2°(k,n) 



4i7T 2 



s n+l 



(13) 



One straightforwardly obtains the following algebraic system of linear equations: 
Q i (x,s)-^(l + e i *)Q°(x 1 s) = a 
- S -(2 + e-™)Q l (x, s) + (l- S -)Q°(x, s) = l-a 



(14) 



which determines the function Q l, °(x, s): 

a^°{x) + b%°{x)s 



Q h °(x,s) 



s 2 + 
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— s — 



a^jx) + b i f{x)s 
p(x,s) 



(15) 



3 + e~ lx + 2e lx 3 + e~ lx + 2e l 



where 



a° a (x) 



9a 



3 + e~ ix + 2e ix 

9(1 -a) 
3 + e~ ix + 2e ix 
3(1 - 2a + (1 - a)e* 
3 + e~ ix + 2e ix 
3a(2 + e~ lx ) 



(16) 



3 + e~ tx + 2e lx 

Denote by s±(x) the roots of p(x, s). Using (|i~3|) one can rewrite 



; /: x 



- ■ 1 r 71 e 
Pz'°(k,n) = — / dx 

Z7T J -TV S i — S 



1 



(17) 



We are interested in the n, k 3> 1 regime, and therefore consider the integrand in ([T7]) for 
x —>■ 0. Here we expose the second order computation, keeping in mind that any order can 
be reached the same way. With 

3 3ix 51x 2 



2 + 20 250 
ix 209x 2 



+ 0{x 6 ) 



(18) 



1 



15 2250 



one gets 
P 3 (k,n) 



2tt 



dx e 



-ikx i 



ix 209x 2 , 



15 2250 



+ 0(x 6 



1 r°° 
2n J-oo 

A 



dx exp 



-ikx — n 



107x 2 ix 



1125 15 



n 



1/2 



cxp 



H25(fc-f|) 2 
A28n 



(19) 
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where A is the normalization constant. 



The expression ([19|) allows one to compute the limiting value of the normalized drift I3 



J 3 = i im Ml 

on the backbone graph TT 3 , where 

/CO _ 
kP 3 (k,n)dk = — 
-00 10 

and, hence, the drift l 3 

i 3 = lim Mi = 2 j = A 

n->oo n 15 

on the graph #3. 

Let us generalize these computations to the case of H q . One can write 

[f]+i 



P fl (*,n) = £ ^(M (2°) 



i=l 



and define the constants pi, 1 < i < [| ] +1 (assuming the existence of corresponding limits): 

Pi = lim lim - g ; ' (21) 



P,(A;,n) 



which satisfy the normalization condition 



[f]+i 

E Pi = 1 ( 22 ) 



i=l 




root 



FIG. 3. Different types of vertices and their corresponding weights (here q = 6). pi gives the 
asymptotic (k, n — > 00) probability of being at a vertex of type i. 
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The sum fl20|) runs over [§] + l non-equivalent vertices (the graph is locally Z 2 symmetric, 
see Fig.^) of the elementary g-gone of the graph. Proceeding in the standard way, we define 
the transform 



Qi(x,s)=T 



P.', 



(23) 



and derive the master equation, whose solution can be expressed in the following form 



P l = T 



— 1 



where the ctj parametrize the initial conditions 



[§]+i 



a,- 



i=l 



and 



/ - 1 

|(2 + T ^-e-) 



M q (x,s,p) 



V 



-pi 

S P3 r —ix 
3 1-/01 e 

s Pi p—ix 

3 1-/91 



8 ga r ~ lx 
3 1 -pi 



§(l + e t!E ) §e' x 

-1 - 

1 3 

- -1 

3 1 






3 





-!+§/ 



For n, k ^> 1 one obtains, using the same method as for q = 3 

P;(k,n)=p l (p)5(k-Up)n) 

where 

= hm <®1 = 

n-»oo n dx 

and s!L is the root of the polynomial det (M g (x, s, p)) the closest to zero. 
To make the system of equations (f24])-(|2"6|) self-consistent we must set 

Pi = Pi(p) 



(24) 



(25) 



(26) 



(27) 



(28) 



(29) 



which closes a system of equations determining pi. Finally, one can write the limiting drift 
in the following form 



L = lim = l q (p) (l -t 



n 



I-Pi 



(30) 



One can check that this formalism gives for q = 3 the same results as has been derived 
above. 



10 



B. Random walk on B 3 : drift and return probability 



1. Analytic results 



We now focus on the braid group B 3 and in particular explain why some statistical 
characteristics of random processes on B 3 have the same asymptotic behavior as the ones 
on PSL(2, Z). The key point is that B 3 is a central extension of PSL(2, Z). Let us recall 
that the center Z of B 3 , generated by a 2 = b 3 = A 2 , is isomorphic to Z. We denote by tc 
the canonical quotient map 



tt: B 3 —> ~ PSL(2,\ 



(31) 



One has then 



7T(o-i; 
7r(cr 2 ) 
7r(a) 



= <J\ 
= o 2 
a 2 



a 2 b 3 
h<x 2 



(32) 



where a 2 and 63 are defined in ([!]) 

A natural representation of the Cayley graph of -B3 is three dimensional. As shown in 
Fig.0|, the map ir can then be viewed as a projection from 3D to 2D. 
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FIG. 4. .B3 Cayley graph and its projection (PSL(2, Z)). Thin arrows correspond to a, thick 
ones to b. Note that 7r(aj) = a, vr(/3j) = /3 and so on. Recall that a>2 has to be identified with a^ 1 

Consider now an n-letter random word w n written in terms of generators of the group 

n 

W n = n a n 

1=1 

where we set cr^ 1 = a_ ri and indices rj are uniformly distributed in {—2, —1, 1, 2}. We recall 
that L B3 (w n ) is the irreducible length of w n . It is evident that 

L B *{w n )>L PSL ^ z \n{w n )) (33) 

(keeping in mind the geometrical interpretation of -B3 shown in Fig.|], we can easily derive 
eq. (|3~3"D following from a triangular inequality). Consider now the irreducible decomposition 
in PSL(2,Z): 

L PSL(2Z) (n{wn)) 

ffK) = II d A (34) 

i=l 
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The asymptotic value of {L PSL<y2 ' z \Tt{w n ))j for n ^> 1 is computed in appendix A by a 
straightforward adaptation of the method introduced in a previous section. From (|3~4]) and 
the definition of the quotient map, we get: 

L PSL(2Z) {n{wn)) 

w n = A 2 ^ n <V 4 ( 35 ) 

where A 2 is 6-letter word in the alphabet {<Ti, cr^ 1 } (see Eq.[|), what implies the following 
condition on f(n): 

71 

\f{n)\ < g 

Hence, the irreducible length of the word w n can be estimated from above 

L B ' A {w n ) < 6|/(n)| + L PSL ^ z \n(w n )) (36) 

Let us show now that the Markovian process f(n) is such that 

(\f(n)\) = 0(y/n) (37) 

1. The symmetry of the process implies that the words A 2 and A -2 appear with same 
probability, which gives (f(n)) = 0. 

2. The increment \f(n + 1) — f(n)\ is bounded from above by some constant. 



Thus, the central limit theorem gives fl37f). This, together with ([33]) and ( p6|) allows one 
to write 



n n n \v n , 

Using the result lim^oo L PSL ( 2 ' z \w(w n ))/n = 1/4 obtained in Appendix for symmetric 
random walk on PSL(2, Z) one arrives at the following asymptotic expression 

/B , = lim ^K) = i (39) 

J n— >oo 77-4 



Statistics of loops on B3: return probability for "magnetic" random walks 

The investigation carried out above shows that if a random walk on the group B 3 ends in 
Z(B$) (we will say in this case a Z-walk), it can be regarded as a closed "magnetic" random 
walk on PSL(2, Z). Namely, if one inserts in each elementary cell of the hyperbolic lattice 
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a "magnetic flux" h (see Fig.fD and denotes by $ the total flux through a closed path on 
PSL(2, Z), then any word w% corresponding to a Z-walk on B3 can be written as 




In other words, the group B 3 being the central extension of PSX(2,Z), gives rise to a 
fibre bundle above PSL(2, Z) such that every full turn around the elementary cell leads 
to another sheet of the Riemann surface of PSX(2,Z). The outcome of this construction 
is that Z-walks on B3 can be decomposed into a product of elementary full turns around 
cells (this is due to the tree structure of the backbone). Hence the function giving 
the probability that a closed n-step loop on a graph PSL(2, Z) carries a flux $ is of great 
interest, especially because at $ = it defines the probability to get a trivial braid (i.e. 
completely reducible word) from a random braid of the record length n. 

First of all we compute m^($) for a walk with local passages in a basis {a 2 , b 3 , a 2 , b^ 1 } 
(let us stress that for magnetic walks a 2 7^ a 2 l ). Denote by §a 2 , ^a^ 1 , $63, jJ&J 1 the total 
number of steps a 2 , a 2 , 63, b% 1 respectively in a given closed path on PSL(2, Z). The fulx 
$ can be written as follows: 



Recall that we consider an n-step process on PSL(2, Z), conditioned by the fact that the 
path is closed (i.e. returns to the origin). Following fl40|) , we rise a simultaneous process $i 
(with $0 = 0) such that 



with 4>i = ±h/2 if the corresponding step on PSL(2, Z) is af 1 , or 0j = ±/i/3 if the step is 
63 x . Evidently the final value $ n gives the total flux $ through the closed path. 

We show that the process $j is not affected by the condition that the path is closed. 

1. Notice that on PSL{2, Z) we have a^ 1 = a 2 , and therefore p(h/2) = p(—h/2). 

2. The sign of the magnetic field can be arbitrarily changed, hence p(h/3) = p(—h/3) 
(i.e. positive, 63, and negative, b^ 3 , elementary turns are equidistributed for closed as 
well as for open paths). 

3. The closure condition on PSL(2, Z) affects the irreducible length of words; the irre- 
ducible forms on PSL(2, Z) being exactly the words of the form a ±1 6 ±1 a ±1 6 ±1 a ±1 ..., 
setting the irreducible length of a word does not change the relative weight of a ±x and 
b ±l in this word. One finally obtains p(h/2) = p(h/3). 




(40) 



i+l — 



(41) 
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The process $j is then a classical one dimensional random walk, and therefore for n large 
one has 

' ,42) 



13 
72 ' 



where = § (± + §) 

This result seems to be interesting in the context of lattice random walks in a transversal 
magnetic field which has relations to the Harper-Hofstadter problem (see [|1^] for review) in 
hyperbolic geometry. 

Returning to the random walk on the braid group in the standard framing 
{<7i, <7 2 , oT , o"^ 1 }, we can compute the distribution for a random process on 

{<7i, <7 2 , erf , a^ 1 }. Modifying slightly the derivation carried out above, one obtains that 
the corresponding process fa is still not affected by the condition of return, and is such that 
p(fa = h/6) = p(fa = —h/6) = 1/2. This yields 

<(»)= ' cxpf-« (43) 

with o-^ = _L. 

The decomposition introduced above allows to compute the return probability, i.e the 
probability p(w n = Id) = p r {n) to obtain a "trivial" braid after n random elementary moves. 
Using ( p5| ) the condition w n = Id is equivalent to the conditions 

L PSL ^ z \rr(w n )) =0&/(n) = 

Denote 

p{L PSL(2,Z) {7T{Wn))=0}=p « {n) 

and 

p{/(n) = knowing L P5L(2 ' z) (vr(w n )) = 0} = ^(n) 
The probabilities Pr{n) and pj;( n ) are independent which allows us to set 

Pr(n) = ^(n)p^(n) (44) 

where p*{ri) is computed in appendix A and p c r (n) can be reexpressed the following way 

p c r (n) = hu°(0) (45) 

Collecting (f4^-(|45|) we arrive at the final expression for p r (n) 

Pr{n) = — (46) 

where A and C are given in appendix. 
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3. Numerical results 



So far there is no constructive algorithm to find the reduced form of words of B3 for 
generators a^. The existence of an algorithm depends crucialy on the set of generators we 
choose. Indeed, it is shown in || that computing the length in terms of generators <7j of 
a braid in B n is an NP-complete problem. Let us mention nevertheless that braid groups 
are "biautomatic" (see 0]) which basically means that there exists a set of generators, for 
which the reduced words are exactly known. This allows in particular to solve the word 
enumeration problem, and to implement methods which can compare two different braids 
in a polynomial time (see [0]). In our case of the simplest nontrivial group B% we tried a 
random reduction procedure, but it converges only in exponential time. Since our analytical 
results are obtained in the regime (n ^> 1), the numerical simulations give no additionnal 
information. 

IV. DIFFUSION ON RIEMANN SURFACES: TRACES AND LYAPUNOV 

EXPONENTS 

We consider the representation of dimension 2 of the groups introduced above, and 
investigate their action on the hyperbolic Poincare plane 7i = {z, Imz > 0}. Namely, we 
consider the following fractional-linear transforms 



We recall that PSL(2, IR) is a subgroup of the group of isometries of TL. The groups 
PSL(2, Z), PSL(2, Z) u , H q , F n admit representations as subgroups of PSL(2, IR) and their 
Cayley graphs (considered in previous section) are now viewed as isometric lattices embed- 
ded into 7i. Now one can investigate their metric properties using the natural hyperbolic 
(geodesic) distance in 7i. We define the lattices under consideration the same way as we 
have defined the Cayley graphs: 

• We construct the set of all possible orbits of a given root point (we choose the point 
i = (0, 1) for conveniency) under the action of the group. 

• We denote by d(w n ) the hyperbolic distance d(i,w n (i)) between i and w n {i). 

• We call "lattices" the Cayley graphs of the groups involved here because of two im- 
portant features: 




(47) 
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— they are discrete subgroups of PSL(2, IR), the group of motion of the hyper- 
bolic 2-space. Hyperbolic distance is a pair-point invariant, that is d(i,w(i)) = 
d(ji, wy(i)), what jusities the term isometric; 

— they have the property of so-called lattice groups: they have no points of accu- 
mulation (for the topology of TC). Recall that for H q , q 6 Z. 



Let us add that the above description is based on well known results on Fuchsian groups 
theory (see p!5|fl^| ). Properties of a Fuchsian group G depend strongly on the fundamental 
domain of G, which is a minimal set of points generating 7i under action of G. The groups 
studied throughout this paper are all Fuchsian groups. We first remind that the fundamental 
domain of the Hecke group is the circular triangle with angles jo, |, - j (see Fig.|5] for H 3 ). 



fundamental domain 
of the modular group 




FIG. 5. Fundamental domain of PSL(2, Z) 



It can be shown that the fundamental domain of F n is a zero-angled n-gone. Our contri- 
bution to this subject concerns the construction of the fundamental domain of the deformed 
group PSL(2,Z) U (Fig.0). 
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T(r) 




FIG. 6. Fundamental domain of PSL(2, Z) u , for u = 1.2. 



We omit technical details of this construction, which can be found in ||14fl . The outline 
is as follows. We first find the fixed point i/u of S u , and xq = 1/(1 — u 2 ) of T u . We 
then draw the only geodesic through i/u which intersects its images by T u and T~ x with 
angle 7r/3. Circles of center xq passing by these intersections complete the construction. 
First notice that the topology of the Cayley graph obtained this way does not depend on 
u. Recall that only commutation relations, independent of u, set the topological structure 
of the Cayley graph). Only the metric properties are affected by u. In particular the area 
of the fundamental domain is finite only for u — 1. The group is then said to be of type 
I in the classification of Fuchsians groups. For u ^ 1 it is of type II. It means that the 
corresponding monodromy problems are deeply different (see |l6j). Solving the monodromy 



problem is an important issue since it allows to get the conformal transform that maps the 
fundamental domain onto 7i. To our knowledge the problem is solved only for u = 1. We 
have therefore to content ourselves with an existence theorem in the general case. Existence 
of such a transform allows to define a map f u from the fundamental domain of PSL(2, Z) u 
to the fundamental domain of PSL(2, Z). The action of PSL(2, Z) u on 7i is in this sense 
conjugate to the action of PSL(2, Z): 

\/lu u E PSL(2, Z)„, lo u {z) = f- 1 o uj u=1 o f u (z) (48) 
The dependence on the parameter u is this way clearly expressed. 



A. Analytic results 

Let us return to the definition of the model and recall that the groups under consideration 
act in the hyperbolic Poincare upper half-plane 7i = {z e C, Im(,2) > 0} by fractional-linear 
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transforms^. The matrix representation of the generators (denoted by hi, 1 < i < n g ) of the 
different groups has been given in section |I[ 

Choosing the point (x ,y ) = (0, i) as the tree root — see fig.§, we associate any vertex 

n 

on the lattice with an element w n = h ak where 1 < < n g and w n is parametrized by 

k=l 

its complex coordinates z n = w n (i) in the hyperbolic plane. 

Strictly speaking Ti should be identified with SL{2, \R)/SO(2); we here identify an ele- 
ment with its class of equivalence of SO (2). The following identity holds (see [p~4| . p~8|] ) 

2 cosh (d(w n )) = Ti(w n wl) (49) 

where dagger denotes transposition. 

We are interested in the distribution function P n (d), and therefore have to look for the 
distribution of traces of matrices w n . The method described hereafter involves mainly the 



results of the paper |13| . The outline of our approach is as follows. We study the behavior of 
the random matrix w n , generated by a Markov chain (which must fulfill ergodicity properties) 
defined as follows: 

w n+1 = M n h an+1 , {a n+1 =i, 1 < i < n g ) with probability — (50) 

n g 

We use the standard methods of random matrices and consider the entries of the 2x2- 
matrix 4-vector V n . The transformation u>„ + i — w n h a reads now 

V„ +1 = I f °, I V. (51, 

This block-diagonal form allows to study one of two 2-vectors composing V n , say v n . 
Parametrizing v n = (g n cos 9 n , g n sin 9 n ) and using the relation d{w n ) = d n ~ 21n^ n valid 
for n ^> 1, one gets a recursion relation v n+ \ = b) a v n in terms of hyperbolic distance d n : 

d n +i = d n + In p a (cos 6) (52) 

where p a is a second order polynomial depending on the specific form of transition matrices 
h a . While for the angles one gets straightforwardly 



2 It is convenient first to define the representation in the Poincare upper half-plane and then use 
the conformal transform to the unit disc. 
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cot 9 n+ \ = h a (cot r , 



(53) 



The action of h n is fractional-linear. 



One has now to study the invariant measure fi(0), giving the asymptotic probability to 
have n = 0. Introducing x = cot 0, we are led to study the action of the group restricted 
on the real line parametrized by x. The statistical properties of /i have been discussed by 
Gutzwiller and Mandelbrot U in the case of the free group A. An alternative, put forward 
in rr ' ' 



13] , is to define /i(x) as the limit of the following recursion relation: 

dh a (x) 



/i (n+1) (x) 



n 



J2^(h a (x 



9 a=l 



dx 



(54) 



The convergence /i^ n '(x) — > /i(x) for n — > oo is assured by ergodic properties of the 
functional transform ( [54] ) and has been successfully checked numerically by comparing to 
direct sampling of different orbits. Despite the absence of rigorous proof, we claim that fi 



is defined with no ambiguity by fl54f) . This enables us to compute the desired distribution 
P n (d). The crucial point required for convergence of /i*" - - 1 to the invariant distribution, is 
the existence of ergodic properties of n . It means that for n ^> 1, the distribution of n is 
exactly given by n(0), independently of n and initial conditions. We introduce the generating 
function for (|52|) ; due to the Markovian structure of (0), we can perform the averaging: 



.. / A (.1 n 



Thus we obtain 



n 



E 

9 a=l 



7r/2 
-7T/2 



[p a (cos( 



lifc 



dOfi(0)[ Pa (cosO)} 



ik 



(55) 



(56) 



This form suggests that for n large P n (d) satisfies a central limit theorem. Indeed such a 
theorem exists (see [|H],^(]]) for Markovian processes provided that the phase space is ergodic. 
We are then led to compute only the first two moments (Lyapunov exponents), which gives 
us 



7i 



(d) 1 ™ 9 
lim = — I dd M mp Q (cosfl) 



n n 



9 a=l 



-tt/2 



and 



where 



2 r ((d-(d)) 2 ) 2 
a = km — \-LL-L = 7a - 7 2 

n— +oo n 



n 



9 a=l 



tt/2 
-tt/2 



dOfx(0) \n 2 p a (cosO) 



(57) 



(58) 



(59) 
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B. Numerical results 



We present in this part the numerical and semi-analytical results for the invariant mea- 
sure \i and the Lyapounov exponent 71. Our main goal is to compare the approach developed 
here with the results following from the study of random walks on graphs (see section |T|. 

Let us call the backbone subgroup B{G) of the group G such subgroup of G whose Cayley 
graph is the backbone of the graph of G. It seems to be more instructive to rely on this 
purely geometrical caracterization of B and to avoid a formal formal definition. Let us stress 
that B(G) is a free subgroup of G. One has for example B(F n ) = F n . Consider now the 
representation of F q by q idempotent generators gi, ..,g q with the following homomorphism 

*-{ Fq ^Z\ u ( 6 °) 
\9i — > b q l a 2 b\ 

Due to injectivity of the following decomposition holds 

H q = U b^(F q ) (61) 
i=i 

with 

^(F ? )nW^) = 0forz^j (62) 



what means that the Cayley graph of H q is the disjoint union of q trees ~W q . Thus we set 
B(H q ) = F q . 

The scale factor Sf is the "average" irreducible length of the generators of B(G) in G. 
In other words, L G {w) ~ SfL B ( G \w) for w G B(G) with L G {w) ^> 1. We have studied 
two different Markovian processes for each group G: (i) simple random walks (characterized 
by the Lyapunov exponent 7*) and (ii) so-called directed random walks (that are walks 
excluding two consecutive opposite steps) on the backbone subgroup B{G) (caracterized by 
Lyapunov exponent 7^). 

By construction ^f/sj gives the average hyperbolic length for an elementary step on G. 
We conjecture that s/7i/7i gives the number of steps to the origin (normalized by n) on the 
graph G. Let us point out that this result links together two definitions of the "drift" for 
random walks on the groups G: the drift / is defined on the graph in metric of words while 
7 is defined in terms of hyperbolic distance for an isometric embedding of G into 7i. Thus 
we claim 
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L G {w)) = ^ (lnTr(W)) (63) 

where a word w is identified with its matrix representation. We believe that equation 
53] is worth interest, since it relates properties of a group defined only trough symbolic 



commutation relations, to geometrical properties of a given representation. 

The stochastic average (...) in ( |6~3|) is necessary, to wash out purely geometrical effects 
such as multifractality investigated in [Tj| . (It corresponds to fluctuations of the hyperbolic 
distance for words of same length on the backbone graph). One has to stress that ( |63"D holds 
due to a "global" spherical symmetry (see |26j for a precise definition of this symmetry for 



graphs) of both models; only the "radial" part of the processes is considered, whereas the 
angular dependence is averaged (here again the ergodic properties play the crucial role). 
This has been checked numerically in the continuous case: generators have to be properly 
normalized, such that each elementary step should have the same hyperbolic length, ensuring 
spherical symmetry, else the invariant measure fi fails to converge. 

All results are summarized in table |I|. 



group 


generators 


backbone subgroup, 
scale factor Sf 


numerical 


semi-analytical 


{d)/n: 

graph approach 


^3 


hi,h 2 ,h 3 


F 3 , 1 


0.3334 


0.332 


1/3 


Fa 


hi,h,2,hi ,h 2 


Fa, 1 


0.501 


0.503 


1/2 


H 3 




^3, 2 


0.1334 


0.132 


2/15 = 0.133.. 


PSL(2,Z) 




^3, 1 


0.2501 


0.248 


1/4 



TABLE I. 
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V. DISCUSSION AND PERSPECTIVES 



We have presented during this work different aspects of random walks on a family of 
hyperbolic groups. On one hand we studied the Cayley graphs of these groups, and briefly 
exposed general methods of computing the Green functions for Markovian processes on those 
graphs; in particular we explicitly calculate the drift in different cases. As an application, 
we studied Markovian processes on the braid group B 3 , and explicitly showed that the drift 
for a symmetric random walk on this group tends at n — > oo to the drift of a process on 
the group PSL(2, Z), which is found to be 1/4. This means that a typical random braid of 
record length n can be released on average by n/4 elementary moves. The graph approach 
and the introduction of "magnetic walks" enabled us also to compute explicitely the return 
probability on B 3 , that is the probability to obtain a trivial (completely reducible) braid 
from a random word of record length n. 

On the other hand we took advantage of the fact that the groups H q and PSL(2, Z) u are 
subgroups of PSL(2, IR) and therefore act naturally in the hyperbolic plane 7i. The Cayley 
graphs of these groups are then naturally embedded in 7i. Instead of the usual length in 
metric of word, we could, thanks to this representation, use the metric structure of 7i and 
study the hyperbolic length of random elements of the group. This problem leads to the 
study of products of random matrices. The method described in |13[ allows us to compute 
the probability distribution of the hyperbolic length. Lyapunov exponents are explicitely 
computed in different cases. 

These two approaches are shown to be related by an equation fl6~3|) . This result is a 
strong motivation for investigating further the geometric properties of hyperbolic groups in 
connexion with other topological invariants. As an example we briefly mention the Alexander 
polynomials. 

The Alexander polynomial Vx(t) of a link K represented by a closed braid w n = ]Jj 
of length n is defined as follows 



[I + t + t 2 ) V K {t) = det 



n 




n 




n 




= det 




+ 1 - Tr 















(64) 



where j runs "along the braid" , i.e. labels the number of used generators, the subscript rj G 
{—2, —1, 1,2} marks the set of braid generators (letters), with the prescription di V = a_j 
and I defines the 2 x 2-identity matrix. For long words (n ^> 1), the following asymptotic 
expression holds: 



1/2 



d(w n )/2 



Tr(w n ) ~ (Tr(w„w^)^ 
One then has, with the parameter u = \f—i (recall that <7j depends on u)\ 



(65) 
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(1 - U 2 + U A ) V K (u) = U 2p{Wn) - u P^n) e d{w n )/2 + j (gg) 

with p(w n ) = Jj(+) — ). In this regime the polynomial is therefore expressed only in 
terms of p(w n ) and d{w n ). The quantity p{w n ) is a "poor" invariant, in sense that it takes 
the same value for a large amount of links. In other words p(w n ) is just the length of the 
element w n projected onto Z. Indeed there exists an obvious group homomorphism tti from 
B 3 to Z defined by Tti{af l ) = ±1. All non abelian properties are lost by this invariant. 
The geometric invariant d(w n ), described above, is much stronger. As we have shown, this 
invariant is related directly to the word length in the group PSL(2, Z), which preserves 
the noncommutative structure of S3 (recall that the random word length in B3 has the 
same asymptotics as the word length in PSX(2, Z)). The information is nevertheless not 
redundant, because there is no nontrivial homorphism from PSL(2, Z) to Z (there is no 
finite order element in Z). In particular, under the condition d(w n ) = (Z-walks), p is an 
exact invariant having sense of a winding number. 

The form (^) seems in particular convenient for possible problems of statistics of Alexan- 
der polynomials, since we know the statistics of both p and d. In particular, for a simple 
random walk, the typical Alexander polynomial could be defined as \7 n (u): 

(1 - u 2 + u A )V n {u) = 1 - e n ^ u)/2 (67) 
where 71 (u) is the Lyapunov exponent of the random product of generators <3"j. 

APPENDIX A: DRIFT ON PSL(2, Z) 

The goal of this appendix is to compute the drift of a random walk on PSL(2, Z) in 
terms of generators ov We keep notations of |II| and proceed the same way, noting that 
the process under consideration is no longer a simple random walk, but is described by the 
transitions shown in fig.[F[ 
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root 



FIG. 7. Random walk on P5L(2,Z) in terms of generators U{. Vertices of type i and o are 
shown. 



A direct counting gives the following master equation for k > 2: 

P*(k, n+l) = \ (P*(k + 1, n) + 2P*(k - 1, n) + P£(k - 1, n)) 
Pg(k, n+l) = \ (P 3 °(A; + 1, n) + P*(k + l,n) + 2P°(k - 1, 



with initial conditions of the form 



P^(k, 0) =aS k , 
Pg(k, 0) = (l-a)4,o 



(Al) 



(A2) 



One then straightforwardly obtains the following algebraic linear system: 



Q l (x, s) (l - S -{e~ lx + 2e ix )j - ^e ix Q°(x, s) = a 
—e-^Q^x, s) + (l - S -{e~ ix + 2e ix f) Q°(x, s) = 1 - a 



(A3) 



determining Q l, °(x, s): 



Q^x, s) 



a^°{x) + b^°(x)s 
p(x,s) 



(A4) 



We omit the details irrelevant for the purpose of this Appendix. We denote as s±(x) the 
roots of p(x, s). They obey the equations 
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s + = 2 - ix + 0{x 2 ) 



s_ = l-- + 0(x 2 
4 v 



(A5) 



and one gets finally 



" 4 



n 



(A6) 



One has now to make sure that for any word w n of n letters in the alphabet U{ the 
following relation holds: 



k(w n ) = L{w n ) + 0(1) 



(A7) 



Even if Fig.|7| makes this statement clear, a more rigorous proof is as follows. Consider a 
given word w, with k(w) = ho- Then the following decomposition holds: 



ko 



(A8) 



with e G {0, 1, 2}, Cj G {1, 2}, G {0, 1}. To prove (|A7|) we use the relation 



L JJ a 2 &3 =A; 



(A9) 



and one has finally 



lim 

n— >oo 



L(w r . 



n 



1 

4 



(A10) 



Let us mention that a more direct derivation of this result can be brought in if one considers 
the PSL{2, Z) generators . The structure of the Cayley graph of PSL(2, Z) depends on the 
basis and in the framing S, T, T -1 it has form of the so-called hyperbolic honeycomb lattice 
(see Fig.R). 
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FIG. 8. The honey-comb lattice. 



Define k — the distance on the backbone graph of PSL(2, Z). The partition function 
P n («) to find the random walker at a distance in k steps along the backbone graph from the 
origin after n elementary steps satisfies the master equation: 

P n+1 («) = ^P n (K + 1) + lp n («) + ^P n (K - 1) (All) 

with the following boundary conditions: 

P„ +1 (0) = ^(P n (0)+P n (l)) (A12) 

This is a standard problem whose solution is known, and the condition L(w n ) = is in 
particular equivalent to k — 0, therefore the probability to obtain a trivial word after n 
random steps (denoted (n)) is given by 

P ;{n) = P n (0) = C-^ (A13) 

with _ 

_ 9 + AV2 . . 2V2 + 1 

C = and A = 

7tt 4 
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